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Abstract 

We study the collapse of a self-gravitating and radiating shell of bosonic matter. The 
matter constituting the shell is quantized and the construction is viewed as a semiclassical 
model of possible black hole formation. It is shown that the shell internal degrees of freedom 
are excited by the quantum non-adiabaticity of the collapse and, consequently, on coupling 
them to a massless scalar field, the collapsing matter emits a burst of coherent (thermal) 
radiation. The backreaction on the trajectory is also estimated. 

1 Introduction 

Much effort has been dedicated to the study of the classical dynamics of gravitationally collapsing 
bodies and a fairly comprehensive understanding of the main features of this phenomenon is now 
available (see, e.g., and References therein). However, it is also clear that classical physics 
is not sufficient for a complete description, and for diverse reasons: first of all, the predicted 
point-like singularity which emerges as the final fate of the collapse is quantum mechanically 
unacceptable, just on the basis of naive consideration of the uncertainty principle; secondly, 
as soon as the collapsing body approaches its own gravitational radius, Hawking radiation Q 
switches on and its backreaction should be included properly |9]. Moreover, the above two issues 
are connected, since Hawking's effect violates the positive energy condition which is a basic 
hypothesis of the singularity theorems 0. One therefore expects corrections to the classical 
picture already at the semiclassical level, i.e., in the region where matter is properly evolved by 
quantum equations on a a space-time whose dynamics is still reliably approximated by classical 
equations. 
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The semiclassical limit for various models has been previously investigated [f| ^, ||, || 
by employing a Born-Oppenheimer (BO) decomposition of the corresponding minisuperspace 
wavefunction [10, 11, |j~2^] . In particular, in Ref. Q the case of a (thin) shell of quantized 
scalar matter collapsing in vacuum was investigated and it was shown that the collapse induces 
the production of matter quanta. Further, explicit conditions were found beyond which the 
semiclassical approximation breaks down. Such a model is, however, both unrealistically simple 
and of little physical use, since the absence of any signal from the shell precludes an observer from 
witnessing the process. In order to overcome the latter shortcoming, an effective minisuperspace 
action for radiating shells was derived in Ref. [13] as the first step in the modelling of a more 
realistic case. This was expected to be useful both for the conceptual problem of unambiguously 
quantizing the system and for obtaining more predictive conclusions: in fact, the possibility of 
adding an observer (e.g., in the form of a detector coupled to the outgoing radiation) allows one 
to define physical (and not just formal) observables [14]. 

As mentioned above, one of the most intriguing aspects of collapsed bodies is the outgo- 
ing flux of thermal radiation predicted by Hawking g]. This effect is usually studied in the 
background of a preexisting black hole, thus separating the problem of the collapse from the 
understanding of the onset of the thermal radiation. Such an approach is inspired by Hawking's 
original computation, where any backreaction on the chosen Schwarzschild background is ne- 
glected and the matter collapsed to form the black hole plays no role, and is further supported 
by the smallness of the (renormalized) energy-momentum tensor of the radiation in the vicinity 
of the horizon ]l~5| . In this framework, one can think of a particle-antiparticle pair being gener- 
ated outside the (event |l6| or apparent ||) horizon, with the positive energy particle escaping 
in the form of thermal radiation and the negative energy antiparticle falling inside the horizon 
decreasing the ADM mass (i.e., proper mass plus gravitational energy) of the singularity. It 
is therefore the horizon, a purely geometrical concept, which appears as the key ingredient for 
this process, and much attention is devoted to studying event horizons as if possessing physical 
degrees of freedom of their own [0, m (see also Refs. ||, || for more recent developments). 
However, on recalling that the central singularity is a transmutation of the collapsed body, one 
is alternatively tempted to define the Hawking effect as a transfer of energy from the ADM mass 
of the collapsing matter to the radiation field. 

The latter statement can be put on a firmer ground if one considers the collapse and the 
onset of the Hawking radiation in the same and only picture, which is also the scheme we are 
naturally led to use if we take the point of view of a (distant) static observer. For such an 
observer the collapse would never end (classically) and the final fate of the collapsing body 
necessarily overlaps with the onset of thermal emission. Hence, the static observer does not 
see an intermediate stage at which there is a true black hole, because the infalling flux of 
negative energy annihilates against the (still) collapsing matter, and it is clear that he would 
consider the Hawking radiation as energy lost by the collapsing body when it approaches its 
own gravitational radius. In this scenario, the applicability of the no-hair theorems becomes 
questionable, with geometry playing (at most) a subsidiary role, and one should be able to 
describe the whole process in terms of dynamical quantities, in much the same fashion as is 



the case for non- gravitational interactions |pl| , and possibly recover unitarity [22, 22] . This 
approach was explored long ago in Refs. [24], whose authors were led to deny the very formation 
of black holes (horizons), and has been revived recently from purely kinematical arguments in 
the framework of optical geometry [25]. To summarize the present situation, it still sounds 
fair to say, as in Ref. g, that the issue of whether the horizon forms or not requires a better 
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understanding of the (quantum mechanics of the) gravitational collapse. 

Instead of constructing a general formalism, in the present paper we shall continue our study 
of a specific (quantum mechanical) model of collapsing bosonic matter, to wit, the self-gravitating 
shell (for a purely classical treatment see, e.g., Ref. |2(|). Given its wide flexibility as a building 
block for more complex configurations, we believe that the conclusions we are able to draw are 
sufficiently general to be taken as hints for other situations. One may query our use of a scalar 
field, rather than fermionic matter, for the collapse. Clearly, in our case one will have, besides 
some simplifications, effects associated with the matter's Bose-Einstein nature, for example the 
formation of a condensate. Nonetheless, even if apparently unrealistic, considerable effort has 
been dedicated to the effects of such matter (e.g., in the context of boson stars p7|]). 

Returning to our previous study of a collapsing self gravitating shell, we consider a "macro- 
shell" constructed from a large number of "microshells" , each of which corresponds to the s-wave 
collapse of a scalar particle of typical hadronic mass (m^ ~ 1 GeV) . In Ref. ^ , we discussed how 
such microshells were bound together to form the macroshell and found the confining potential 
to be linear and dependent on the shell radius. This implied that as the macroshell collapses, 
the confined microshells undergo non-adiabatic transitions from the ground to higher excited 
states. The attitude we adopted is that, because of their bosonic nature, the microshells formed 
a "condensate" and that, once enough microshells were excited (thus widening the macroshell) 
they would collectively decay to the ground state by creating additional particles (in analogy 
with hadronic string theory wherein a long string - excited hadron decays into short strings 
- lighter hadrons). Thus during the collapse the creation of additional microshells (particles) 
leads to a backreaction slowing the macroshell. 

It is clear that in the above the Schwarzschild (ADM) mass of the shell does not change. We 
previously indicated that a more realistic approach would be to consider non-adiabatic effects 
leading to transitions from the ground to the higher (excited) states for some of the bound 
microshells and, subsequently, these states would decay and the microshells would return to the 



ground state by emitting radiation |TJ, 14|. At the end of such a process, the proper energy 



of the macroshell would be essentially unchanged (except for the small change in ground state 
energy due to the change in macroshell radius) and the net balance is then a transfer of the 
gravitational energy of the macroshell to the radiation field which decreases the ADM mass of 
the system and modifies the trajectory for the radius of the macroshell. 

In our present approach we do not consider particle (microshell) creation, which would 
require much greater energy than the excitations of the microshells due to non-adiabatic effects 
(in Ref. ||] we also showed that microshell creation occurs for rather extreme cases, in which the 
size of the horizon is comparable to the Compton wavelength of the microshells), but study the 
latter excitations. In order to determine the latter excited levels we must describe the motion of 
one microshell in the mean gravitational field of the others and obtain an effective Schrodinger 
equation for the microshell. Such a study is attempted in the next Section. Subsequently 
one must determine the amplitude for the excitation of a microshell, as a consequence of the 
time variation of the binding potential due to the collapse of the shell (this is attempted in 
Section |3.1[ ). The model is constructed bearing in mind a number of constraints: firstly, since 
our scalar particles have a Compton wavelength of ~ 10~ 14 cm, one must have both that the 
radius of the shell and the thickness of the macroshell be much greater, otherwise fluctuations 
due to the quantum nature of matter will destroy our semiclassical description of the macroshell 
motion (see the consistency conditions for the semiclassical approximation in our previous paper, 
Ref. H). This essentially implies that the number of microshells of hadronic mass be at least 
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10 so that, for example, the gravitational radius is sufficiently large. Another constraint which 
will allow us to simplify our results is to suppose that the collapse be non-relativistic (i.e., "slow" 
all the way down to the horizon), which implies that only the first microshell excited states will 
be relevant. 

In the remainder of Section || we address the emission of the radiation (scalar massless 
particles - we shall subsequently just call them photons) and the macroshell trajectory with 
backreaction. Again we shall enforce some constraints, in particular that the coupling of the 
radiation to the microshells be sufficiently large so that the decay time is smaller than the 
unperturbed collapse time. Actually, we shall see that such is already the case for a very small 
coupling constant. Further, interesting features of the model are that the collapse is slow and the 
shell tickness remains much smaller than the typical wavelength of the emitted quanta, so that 
many emissions occur practically in phase, thus rendering the radiation process highly coherent 
and the backreaction great. As a consequence, it will be sufficient for us to just consider first 
order perturbation theory in the radiation coupling constant and keep terms to lowest order in 
the macroshell velocity. 

One may wonder if our model is useful in describing Hawking's radiation, which has two main 
features: the first one is the thermal spectrum and the second one is its "observer dependence" . 
The former follows from the adiabatic hypothesis that the ADM mass of the source (black hole) 
changes slowly in time and, therefore, one does not necessarily expect to recover a thermal 
spectrum in a fully dynamical context in which the backreaction is included. By the latter we 
mean that a freely falling observer remains in its initial state all the time and should thus not 
experience the flux of outgoing energy measured by a static observer, regardless of the adiabatic 
approximation. We just point out that our model shows such a sort of "observer dependence", 
since the proper mass, which remains constant, is the energy as measured by the observer 
comoving with the shell which, in turn, does not experience any change. At the same time there 
is a flux of outgoing radiation associated with the decrease of the ADM (Bondi) mass. However, 
whether such a radiation can be consistently identified with Hawking's or should be considered 
as a totally different effect will need further inspection. Let us note in any case that, in contrast 
with Hawking radiation, the coupling constant of the radiation to the microshells will appear in 
the final expressions. 

Naturally our results will also be constrained by our availability of computing power: bearing 
all the above points (and approximations) in mind, our results are illustrated and summarized 
in the Conclusions. 

We use units for which c = 1 but explicitly show both Newton's constant G and Planck's 
constant h. The Planck length and mass are then given by £ p = TiG and m p = h/t p . 

2 The macroshell inner structure 

In order to study realistically the evolution of a collapsing body it is important not to neglect its 
spatial structure. For instance, the equivalence principle implies that a point-like freely falling 
observer does not experience any gravitational effect. However, any object has a spatial extension 
and tidal forces become effective when this size gets close to the typical scale of variation of the 
space-time curvature. At the same time, it is also important to keep the model simple so as 
to allow one to study it. A good compromise between realism and simplicity is given by the 
macroshell which was introduced in Ref. || and which we now review and further develop. 



4 



We view the collapsing shell not just as a singular spherical surface in space [£8[ but as a set of 
a large number N of microshells which correspond to collapsing s-wave bosons and are described 
by a wavefunction (thus the microshell radius is actually the "average" radius). Each microshell 
has negligible thickness and proper mass m such that their total proper mass equals the proper 
mass of the macroshell, Nm = M. If we denote by r £ (0, +00) the usual "areal" radial 
coordinate, we can order the microshells according to their area, so that r\ < T2 < ... < rjy, 
and assume that the thickness 5 of the macroshell is small, 

< 6 = tat - n < n . (2.1) 

The space between each two microshells has Schwarzschild geometry, 

, 9 ( 2GM,i\ , n / 2GM i \~ 1 l2 2 o . , 

dsf = -(l — - Jdtf + fl — -J dr 2 + r 2 dn 2 , (2.2) 

with a suitable ADM mass Mj < M i+1 , M = for r < t\ and M N = M s (total ADM mass) 
for r > rjv- For this description to be consistent, one must check that the microshells do not 
spread during the collapse and 5 remains small in the sense specified above. This was already 
shown in Ref. || to be ensured by the mutual gravitational interaction of the microshells if N 
is sufficiently large and the average radius of the shell R is much larger then the gravitational 
radius Rh = 2 G M s . We shall provide further evidence for this in the following. 

Before studying the dynamics of the microshells, one must also "fix the gauge" corresponding 
to the freedom of making coordinate transformations on the relevant solution (|2.2j) of the Einstein 
field equations. This can be done by choosing an ADM foliation p9f | for the corresponding space- 
time manifold. Two options which cover the domain of outer communication (the region outside 
the horizon) are given by the hypersurfaces {S r } of constant proper time r, associated with 
observers comoving with the microshells, and the hypersurfaces {S^} of constant Schwarzschild 
time (denoted by too = i/v) which are associated with (distant) static observers []. Since we 
shall consider only bound orbits for the shell with initially large radius and negligible velocity, 
it is useful, in this respect, to separate the collapse in the Newtonian regime (NR) and the near 
horizon regime (NHR): 

NR) when the radius of the shell is large (and its velocity small), the difference between the 
proper time of the microshells and the Schwarzschild time is negligible, meaning that 
observers comoving with microshells at different radii can synchronize their clocks with 
the clock of the distant observer. Therefore, one expects that the Hilbert space which 
contains states of the microshells as viewed by the comoving observers, 7i T = {<J>( T )(r, r)}, 
is the same as Tit = {^\r, too)} viewed by the (distant) static observer. We examine this 
case in Appendix |A|. 

NHR) near the gravitational radius, the relative gravitational red-shift between different parts 
of the macroshell is no longer negligible and the two foliations {E T } and {S^} differ sig- 
nificantly, since too = £00 ( r , T ) is a large transformation there [see Eq. ( |B.4| )] . One can still 
introduce an Hilbert space TL T as relevant for the comoving observers, and a space Ht for 

1 Strictly speaking, the domain of outer communication is well defined only provided one knows the development 
of the whole space-time manifold Q. However, in our case we cannot know from the beginning whether an event 
horizon will form or not and, thus, if {Et} can cover just a portion or the whole of space-time. At the same time, 
{£ r } is defined naturally just on the world-strip where the macroshell has support. 
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Figure 1: The macroshell structure. 

the distant observer, but the relation between the two spaces (or their physical equiva- 
lence) is far from trivial, as will be apparent from the corresponding Schrodinger equations 
we display below and in Appendix || (see, e.g., Refs. |3^, |23|] for analogous considerations). 

In NR there is no conceptual ambiguity, since the Newtonian (Schwarzschild) time (equal to the 
proper time) is a Killing vector with respect to which canonical quantization is straightforward 
and most calculations can then be completed analytically (see, e.g., Appendix ^] and Ref. M). 
Instead, in NHR there is no a priori reason to prefer either of the two Hilbert spaces to describe 
properly the matter in the shell. Further, if one were to neglect the backreaction on the metric 
and study a shell evolving on a fixed (Schwarzschild) background, the Schwarzschild time would 
still be a Killing vector and provide a unique way for canonical quantization. However, if one 
considers a collapsing macroshell with a microshell structure in order to successively include 
backreaction, the canonical analysis is much more involved and ambiguous. 

In the spirit of locality, it seems natural (and also appears much easier) to analyse bound 
states of microshells with the aid of 7i T and then move on to the point of view of the static 
observer by simply making a change of time coordinate (after having taken the thin shell limit). 
However, in so doing one does not expect to recover the space Tit obtained prior to the thin 
shell limit, nor is it obvious whether the Principle of Equivalence, in any of its forms plfl , can be 
trusted for length scales associated with such bound states which are of the order the Compton 
wavelength of (elementary) particles. In the present paper, we shall motivate taking the thin 
shell limit at some point by studying the trajectory of the shell only down to a radius greater 
than the gravitational radius (horizon) by several shell thicknesses. 

In order to define an effective Hamiltonian for each microshell, let us single out one of the 
N microshells, denote its radius by x and study its motion in the background defined by the 
remaining N — 1 microshells (see Fig. ||). When x 7^ r^, one can apply the junction conditions 
|28| which yield 



dx\ 2 (M>-M<) 2 M > + M< G 2 m 2 _ 

dTJ =~ 1 + ^ + G — x — + ^ = F{xh (2 ' 3) 



where M< (M> = M x ) is the ADM mass computed on the inner (outer) surface of the microshell 
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Figure 2: The equivalent configuration when the microshell lies outside the macroshell: 
x < X; (2) X < x. 



and t x is the proper time of the microshell. Then, on multiplying Eq. Q2.3] ) by m/2 one will 
obtain the equation of motion for a microshell in the form of an effective Hamiltonian constraint 

1 / dx \ 2 

H m =-m — +V = . (2.4) 



2 \ dr x / 

Since m <C M, we shall take M> — M < ~ m, which amounts to dx/dr x ~ for large x (S> Rh) 
and is consistent with the choice of having the microshells confined within the thickness 5 at 
large radius. With this assumption both M> and M < become functions of the microshells 
distribution {rj,Mj} only, but the evaluation of V remains extremely involved because one 
should first determine the most likely set of 2 (N — 1) values {r^, Mi}. 

One can get an approximate expression for V by making careful use of the thin shell limit 
(5 -C r%) and considering the form of the potential when x < r\ or x > n + 5. In this case, 
if we denote by X the (average) radius of the macroshell of proper mass M — m, we can also 
supplement Eq. (2.3) with the analogous expression for the macroshell. For x < X [see case (1) 



in Fig. |2| one then has 



dxV 



Gm G 2 m 2 

+ —^ = F < 



dr J x AlX 2 



(2.5) 



dr J \M-mJ X AX 2 <v ; 



where we have equated the proper times of the two shells in agreement with our choice of 
foliating the space-time into spatial slices S T parameterized by the local time of the microshells. 
Analogously, for X < x [see case (2) in Fig. |2| 



dx\ 2 _ G 2 M s -m + G 2 m 2 _ ^ 
dr J x 4x 9 



(2.6) 
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One is now ready to define the relative and centre-mass radius according to 



r = x — X 



R = T7 X + -TP X 

M M 



(2.7) 



where fi = m (M — m) /M is the reduced mass of the system. The above relations can then be 
inverted and, upon substituting into Eqs. (2.5) and (|2.6|), one obtains an effective Hamiltonian 
for the two shells given by 
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dR\- 
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2 \dr J 2 P \dr 
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V (r) = H (r) + H (r) _ 



(2- 



The potential contains the following relevant terms 
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f< "2 ' 



(2.9) 



where we have also assumed m <C M s (< M) and neglected all non-leading terms in m. The 
assumption M s S> 77i, although is not guaranteed a priori, has been done here just for the 



purpose of simplifying the displayed expressions. However, we shall show in Section 3.4 that it 
actually holds true all along the collapse in all the cases of interest. 

At this point, one can use the external potential to shape the distribution of microshells 
inside the macroshell and construct the potential for |f| < 5/2 (see Appendix [A| for an 
explicit calculation in NR). By expanding the latter potential in powers of 5/R and neglecting 
terms of order (5/R) 2 or higher, one finally obtains an effective Hamiltonian for each microshell, 



H (r) 



1 



m 



dr 

7h 



(2.10) 



with a potential function which interpolates smoothly between the two linear in r (outer) parts 
of the potential given by Vi r) in Eq. 



y rn. 



GM s m 
2R 2 




G 2 M 2 m 
AR 3 



\r\ < 



(2.11) 



We then observe that for R 3> Rh ~ 2GM the term linear in f is negligible and the potential 
is thus symmetric around the average radius of the macroshell, but for M s < M it becomes 
steeper for f < and flattens out for f > (see Fig. || for the case of N = 10 40 protons and 5 of 
the order of the Compton wavelength £h of the proton). The same effect happens for decreasing 
values of the average radius at fixed ratio M/M s (see Figs. |] and ||; note that, for this choice of 
M, the potential confines all the way down to Rh)- 



Figure 3: The potential for M = 10 40 m h , R = 4R H = l(T 13 m, <5 ~ l h = 1(T 16 m, 
— 5 < f < 5 and 1 < M/M s < 5. Vertical units are arbitrary. 




Figure 4: The potential for M = M s = 10 40 m h , 5 ~ £ h = l(r 16 m, -5 < r < 5 and 
Rh < R < 2Rh. Vertical units are arbitrary. 
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Figure 5: The potential F (t) for M = M s = 10 40 m h , 6 ~ £ h = l(T 16 m, -5 < f < 5 and 
-Rf//4 < R < Rh/2- Vertical units are arbitrary. 



For a comparison with the case when the space-time is foliated into slices parameterized 
by the Schwarzschild time measured by a (distant) static observer we refer to Appendix [B| 
There we show how to obtain the relevant potential Vm as an expansion in powers of G. Then, 
since Vm contains terms of all orders in G, while V^P goes only up to 0(G 2 ), this comparison 
clearly confirms that perturbative methods cannot always be trusted when dealing with strong 
gravitational fields or, if we wish, the disadvantage of using the Schwarzschild rather than proper 

(t) 

time. In the same Appendix we show that, if one takes the this shell limit of V m and then 
changes to the Schwarzschild time multiplying by [dr / dt oo) 2 , one recovers the thin shell limit of 
V$ in NR. This means the two steps of introducing a potential and taking the thin shell limit 
commute (as expected) in NR. 

For R S> Rh and M ~ M s , the potential V m is symmetric around r = 0. This confines 
the microshells within the following thickness ||] (see also Appendix |C] for more details on the 
spectrum of bound states) 



£m. I R 



4 /3 W#) , (2-12) 



(2JV)V3 \£ p J m \R H 

where £ m = h/m is the Compton wavelength of a microshell. For R close to Rh or smaller ratios 
of M s /M the potential twists (see Figs. || and[3|) and, eventually, ceases to confine the microshells 
on the outside if the macroshell shrinks below the gravitational radius [see Eq. ( C . 1 4j ) ] . We 



further note that the spreads of the lowest wavefunctions are also of order 5, thus the bosonic 



2 We assume there is no discontinuity experienced by the test microshell when it crosses the borders of the 
macroshell and equate the first derivative of the quadratic potential to the first derivative of the linear potential 
at \f\ = 8/2. 



10 



microshells are essentially superimposed and form a condensate, hence some classical singular 
behaviour, such as microshell crossings, do not occur. 

Let us now summarize what we have done in simple quantum mechamical terms. Since 
we are considering wavefunctions, we have viewed each microshell as immersed in a continuous 
matter distribution, corresponding to the mean field of the others which are then treated as a 
(wide) single shell (Hartree approximation or independent particle model - see also the Thomas- 
Fermi model of the atom |32| , ||). Then we noted that, inside the macroshell, a microshell can 
undergo small displacements (S/R si small) about equilibrium, which will be harmonic, and 
on requiring continuity with the potential outside the shell one obtains Eq. ( |2.11| ). Further, 
since in the context of our Hartree approximation the total wavefunction is just the product of 
the individual ones f], it will be sufficient in what follows to just exhibit the Lagrangian and 
Hamiltonian for a given microshell. 



3 The radiation and backreaction 

In this Section we use the description of the macroshell as given in Section [2] and add the coupling 
to an external radiation field. This will allow the shell to radiate away any excess of the proper 
energy induced by the non-adiabaticity of the collapse, as outlined in the Introduction. 

Our aim is to evaluate the flux of emitted radiation, which will determine M s = M s (r), 
and the corresponding backreaction on the trajectory R = R(r). We shall first estimate the 
non-adiabatic amplitude of excitation and, subsequently, give a purely quantum mechanical 
(coherent) treatment of the entire process of excitation and emission by coupling explicitly the 
microshell modes <3? n to the radiation field. The outcome will be a set of two coupled ordinary 
differential equations (for R and M s ) which we shall solve numerically. 

3.1 The excitation amplitude 

From the previous Section we know that the microshells are governed by the Schrodinger equa- 
tion (see also Appendix [C]) 



~d~T 
with 

-Jl 



i^ = ^$, (3.1) 



H$ = ^~ + VW , (3.2) 

2 TO 



where Vm is given in Eq. ( 2.11 ) and has an explicit time dependence due to R = R(t). This 



type of equation can be solved by making use of invariant operators / = J(r) which satisfy |p3| 

ai 



. (3.3) 



The general solution <£ s = $ S ( T ) can then be written in the form 

\^,r) Is = Y,c n e l ^\n,r) I , (3.4) 



3 Let us also remember that the many-boson wavefunction is symmetric. 
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where | n,r) I is an eigenvector of / with time-independent eigenvalue X n and the c n are complex 
coefficients. We also recall that the phase ip n = (p n ( T ) is given by the sum of the geometrical 
and dynamical phases, 



= t I i {n,T' \Ud T > +iH$ \n,T l ) I dT l 



(3.5) 



The Hamiltonian in Eq. ([O]) describes a harmonic oscillator of fixed mass m and variable 
frequency 



o 1 Rh 

Hence, one can introduce standard annihilation and creation operators a and a\ 



(3.6) 



~2h 



f + 



m £1 



with 



a, a 1 



1 and define the vacuum state as 

a|0,r) a = 

A complete set of eigenstates A = {\ n,T) a } is then given by 

(».)■ 



0,r> a , 



(3.7) 



(3.8) 



(3.9) 



where 



a\ n,r) a = yfti\ n-l,r) a 



(3.10) 



a' n 



,r) a = Vn + l| n+ l,r) a 



The solutions of the corresponding time-independent problem obtained by assuming R constant 
and displayed in Appendix |C] are here recovered as <&n\r) = (f | n, 0) Q . 

One can also introduce the linear (non-Hermitian) annihilation and creation invariants b and 

6t mm, 



'2% 



+ i (X7r r — mxr) 



(3.11) 



where 



k St 



1 and the function x = x(t) is a solution of 



x, + Q x 



m 2 x 3 



(3.12) 



with suitable initial conditions. The system then admits an invariant ground state defined by 



b\ 0,r) 6 = , 



(3.13) 



4 A dot denotes derivative with respect to r throughout the paper. 
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from which one can build a basis of invariant eigenstates B = {| n, t) 6 }, 

\n,r) bs ee 

where 



n, t 



> ' /6 > 



(3.14) 



S| n,r} 6 = v^l n,- l,r) 6 
6' | n, t ) 6 = \Ai + 1 1 n + 1, r 



(3.15) 



6 ' 



We can then introduce the invariant number operator (in analogy with the standard number 
operator fit fi) 



The two basis A and B are related by Bogoliubov coefficients according to 
fi = B* b + A* St 



where 



it = Btf +A"b 



1 



6 = 5 fi - A* fit 
St = B* fit - A a 



(3.16) 



(3.17) 



A = - VmQ 



1 



B = - VmU 
2 



1 

m SI x 



n 



X 



(3.18) 



Further, A and ,B will coincide at r = if 



6(0) = fi(0) 



nn ( j c (o) + -j = # m (o) 



(3.19) 



that is -A(O) = and B(0) = 1, which means that for the function x we must require 

1 

x(0) 



v/mfi(0) 
x(0) = . 



(3.20) 



Since $7 ~ i? and -R(O) ~ 0, an approximate solution to Eq. ( [3.12 ) compatible with the initial 
conditions in Eq. ( |3.20| ) is given by 







2Vm ft 3 



(3.21) 



provided r is sufficiently short so that x is negligible. The corresponding Bogoliubov coefficients 
are then given by 



A ~ 



(3.22) 
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and satisfy the unitarity condition 



\B\ 2 - \A\ 2 = 1 . (3.23) 

The amplitude of the transition from the state $0 to after the time r > can be easily 
computed. In particular, one finds []33| 



l ^2n(T) = a (2n,r|0,r) b 



1 MTVW (3 . 24) 



Vfi V -B J 2 n n! 

A)-2n+i(T) = a (2n + l,r I 0,r) 6 

= . (3.25) 

To leading order in \A\ (i.e., $7 ~ R), one then has Aq^q(t) ~ 1 and the amplitude of the 
transition to the excited state with energy E% n = Eq + 2 n fi 0, is 



i 2 (r) ~ t nVWW^(rl 



= H) n ^#rf^) (3.26) 

V > 3 n 2 n/2 n | W ' V > 

where we have assumed Rh = 2GM S remains constant in the interval (0,r). The relevant 
expression for us is thus given by Eq. ( p. 26 ), in which we recall that R < 0. We also remark that 



n/2 



the above transition amplitude follows as both a non-adiabatic effect (since A^2n oc R n ) and 
a consequence of the finite thickness of the macroshell (since Aq^u 5 n ^ 2 ), the latter being 
further related to the quantum mechanical nature of the bound states [since 5 oc 

An order of magnitude estimate for the most probable transition (2 n = 2) can be obtained 
by setting R 2 ~ G M s /R and using Eq. (2.12), in which case 



(I V /3 1 

\A ^ 2 (R, N)\ 2 < \A ^ 2 (R H )\ 2 ~ - ^ , (3.27) 

where we also approximated M s as N m in order to obtain a function of just the number iV 
of microshells. Further, one can also check that the above probability is less than one and, for 
realistic cases, quite small (we shall exhibit this later). This implies that the probability for 
a microshell to get excited a second time (after it has radiated away the energy E 2 once, see 
next Section) is negligible and, considering also that the emission probability is small, we shall 
not consider double emission. We further note that, as expected, the larger Rh the smaller the 
"tidal" effects and the excitation amplitude. 



3.2 The (thermal) radiation 

The problem of coupling a radiation field in the form of a conformal scalar field to collapsing 
matter was previously treated in Ref. || (see also Ref. 0) for the case of a collapsing sphere 
of dust, for which we recovered the Hawking temperature in a suitable approximation. For 
the present situation, we shall find it more convenient to use the analogy with an accelerated 
observer in flat space-time as put forward in Ref. j3f| • 
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The basic observation is that, if the proper mass or the ADM mass (or both) vary with time, 
the collapse of the shell is not a free fall. This can be seen very easily by noting that the solutions 
to the equation of motion Q2.3]) when M > and/or m depend on time differ from the case when 



the same parameters are constant. Alternatively, one can use the second junction condition [28 
to compute the surface tension V of the shell. As we have shown in the previous Section, the 
proper mass M of the shell increases in time, therefore, although our choice of initial conditions 
for R, M and M s are such that 'P(O) = 0, the tension subsequently increases when M increases. 
The tension then slows down the collapse S, and this effect should be further sustained by 



the emission [13], which very quickly brings the proper mass back to the initial value. As we 
outlined in the Introduction, we shall then consider M as effectively constant along the collapse 
and compute the net variation of M s in time. 
Let us denote by 



R-R 



free 



(3.28) 



the difference between the actual acceleration R of the shell and that of a freely falling body, 
Rfree, having the same radial position R and proper velocity R. Of course, a would equal the 
surface gravity, 



1 



a H 



AG M H 



(3.29) 



of a black hole of mass M s for R = Rh constant. In this limiting case, one could exploit 
the analogy between Rindler coordinates for an accelerated observer in flat space-time and 



Schwarzschild coordinates for static observers in a black hole background [15] and find that the 
shell emits the excess energy with the Hawking temperature 



h an 



(3.30) 



where ks is the Boltzmann constant (an explicit construction which leads to this result was 
proposed in Ref. ||). In general, however, the effective acceleration a will have no a priori fixed 
relation with M s and one should keep it as an independent function of r to be determined later 
consistently with the equation of motion of the shell. It is also clear that, should R turn out to 
be constant, a will correspond to the acceleration associated with the instantaneous ADM mass 
and distance from the horizon. 



We then consider an isotropic massless scalar field tp conformally coupled to gravity [15] and 
to the microshells. The Lagrangian density for a given microshell and the radiation field is given 
by 



l 



+12 E 



n>0 



2 y^n^n ^n^nj 2 m dr dr ^n^n ^ n !■ ^ 

1 



(3.31) 



where 1Z is the curvature scalar and Vm^ is the potential given by Eqs. Q and The 



last term in the first square bracket represents minus the interacting Hamiltonian Hi nt with e 
the coupling constant of the radiation to a microshell while the wavefunctions of the latter are 
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expressed in terms of the states found in Appendix^] as $ n /r, where we exhibit the energy levels 
{n > 0) so that the form of the interaction can be chosen in such a way that - to lowest order - 
it just contributes to the transition between different levels and does not modify the macroshell 
ground state energy. Further, the factor of (1/r) comes from the normalization measure now 
being r 2 dr. In fact, in the thin shell limit it is consistent to approximate the metric on the shell 
with the outer Vaidya line element |3(| [37]]. Let us again note that we have just exibited the 
Lagrangian for a given microshell since in the Hartree approximation the total wavefunction is 
just the product of the individual ones 5 . 

Another important observation comes into play at this point. In the preceding Section we 
have computed the (non-vanishing) transition amplitudes ^4o-*2n to lowest order in R, that is 
O(R) [see Eq. ( 3.26| )1, with M and M s held constant. Therefore, in order to determine the flux 



of radiation to lowest order in R, it is sufficient to compute the probability of emission to O(R ) 
with M and M s constant. In this approximation (M s constant), the Vaidya metric reduces to 
the Schwarzschild metric, so that the relevant four-dimensional measure of integration is 



^/^cfx~r z sin z 6d6d(j)drdt , (3.32) 

where t is the Schwarzschild time measured by a static observer (denoted by t^ in previous 
Sections) , 

(3.33) 



VI - Rh/R 



Finally, since the functions <£ n are spherically symmetric and peaked near r = R(t), it is 
convenient to integrate over the angular coordinates and write 



y /=gd i x = Air{R + fydfdt, (3.34) 

where f is the relative radial coordinate and R = R(t) is the trajectory of the macroshell as 
before. 

After the lapse of proper time r, each microshell will jump into the excited state of energy 
E2n = EQ + 2nfi£l with a transition amplitude Ao_ >2n( T )- Successively, it can decay back to the 
ground state by emitting a quantum of energy hiu (having a wavelength much larger than the 
shell thickness, see Table [l]) of the scalar field. We shall later see that the macroshell velocity 
is small and the emissions are numerous. As a consequence, the distance covered between each 
emission is small with respect to the wavelength and therefore the emitted radiation will be 
coherent. 

Let us now estimate, using first order perturbation theory in e, the transition amplitude for 
the emission of quanta of the scalar field, which will occur when the microshells in the state <3?2n 
decay back to the ground state $n, f° r the interval between the times t\ = t{j\) and ti = t(r2) 
(0 < r\ < t%). It will be given by 



/ dt ( final | H mt \ initial ) = — — \] / dt ( l w ; t \ <p | 0; t\ ) 1 | (j,0;t 
Jt! n i=1 J tl j=1 



N 

X EE%^EI1 \r,s;t)(r,s;t\r,0;0) 

n>0 l^n s>0 r=l 



J Of course, both the total and the single particle wavefunctions are normalized to unity. 
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47rie 



N 



h £ I ' dt J df i $ o(rO * 2 (fi) A ifi ^ 2 (t, 0) 
x e - 2<nT (l w I^Cfl + fi.t)! 0) , (3.35) 

where (i,n;t \ i,0;0) = A^o— > n (i) is the excitation amplitude (for the i th microshell) which, as 
computed in the previous Section, is dominated by the contribution with n = 2 and for the 
photon state | ) we consider the Unruh vacuum, corresponding to an outgoing flux of radiation 
[15]. This of course implies that a is non-zero, indeed we shall actually find that a ~ — i?f ree (since 
R ~ 0), consistently with our choice of vacuum and, near the horizon, the flux corresponds to 
Hawking radiation (see Ref. 3J] for a detailed discussion of the analogy between Schwarzschild 
and Rindler spaces). Let us further note that, although a black hole is only close to being 
formed, it is generally accepted that it is the Unruh vacuum that best approximates the state 
that would be formed following the collapse of a massive body |15| . 

On taking the modulus squared of Eq. ( |3.35 ), one obtains the transition probability 



16 7T 2 e 2 ^ ^ 

i=l j=\ 



*2 



dt" 



t-2 



dt' I df'l / dr'j 4' j0 _> a 4- _ 2 . 



-i2(n" t"-Q't i ) 



: $ (rf) Mr") Mr'i) Mr') ( | f>(R" + ?( , t") 0(R' + f',t>) | ) 



(3.36) 



where X' = X(t'), X" = X(t") for any function of time and the last term is the Wightman 
function for the scalar field <p. 

The four-dimensional Wightman function can be related to its two-dimensional counterpart 
[15] through 

( | tp(R" + f'i t") <p(R' + f' p t') | ) 



D+(R" + f'/,t";R' + r' p t" 



4n(R" + f?) (R' + f'j) 

H In [(Avij — i e) (Auij — i e)] 
16 vr 2 (R» + ff) {R> + f' ) 



(3.37) 



where, of course, we only consider s- waves and 

Av i: j = 

~ At + 



t " - t ' + {R" + f'(),-{R' + f' j ) 
RjAR + Afjj) 



R-R H 



(3.38) 



In the above r* = r + ln[(r/.R#) — 1] is the turtle coordinate p9|. We also set At = t" — t' 



(similarly for AR and Af), R = (R" + R')/2 and assumed R - R H > |Ai?|, [Af. 



' 31- 



Further 
An 



Kl 



-2Rh |exp 



t" - (R» + r>>) 



2 Rh exp 



2Rh 
T-h\ I 



exp 



f- (P' + f^ 



2«h 



exp 



exp 



2R H 

1 /At gjgj-Ag/2) \ 
2P^ ^ 2 / 
1 / R M + AR/2) At s 



2R H 



R- R H 



(3.39) 



17 



where T = (t'+t")/2 — t\ and we have performed the same approximation as we used in obtaining 
Eq. 

The above expressions can be simplified on noting that one expects At 3> 5, AR (= 
R At y/1 — Rh /R) corresponding to the fact that the wavelength of the emitted scalar quantum 
is much greater than the width of the shell and the velocity of the shell sufficiently small so that 
it is not displaced much between one emission and the next (see Table Q). One then obtains for 
Eqs. ( 3.38 ) and ( 3.3S| ) respectively 

Av ~ At (3.40) 



Au ~ 2Rh exp 
= 4:Rh exp 



R* — T — t\ 
2R^ 

R* — T — t\ 
2R~h 



exp 



sinh 



At 
At 



exp 



At 



(3.41) 



On introducing L = ti — ti, expressing Eq. (|3.36| ) in terms of T and At and approximating 
the integration measure as in Eqs. (3.40) and ( 3.41| ) with 



R' = R(T) - R(T) Jl - R H /R (At/2) 



(3.42) 



R" = R(T) + R(T) sjl - R H /R {At/2) , 
and R(T) = R, one obtains to leading order in R 



P(2Q;ti,L) 



N 



2 e 2 h f L jrr \A ^ 2 \ 



dT- 



2 m 2 Jo R e Q 2 



d(At)e 



x In 



AR H e {R *- T - tl)/2RH (At-ie) sinh 



At 
4i^ 



i e 



dT- 



I A 



0^2 1 



e 2 h 

2 m 2 Jo R 6 n 2 
H) 

20 i/l - R H /R 



d(At)e 



1 1 

+ — — — coth 



At-ie AR H 



At 

Tr~h~ 



l € 



(3.43) 



with A = L — 2\T — L/2\. In the above, we approximated r" — r' ~ At \J\ — Rh / R according to 
Eqs. ( |3.33j ) and the last equality follows from integration by parts in At (neglecting boundary 
terms) . 

The integral ( p. 43 ) can be computed by analytic continuation in the complex At plane in 
which, upon using the formula 



coth(7rx) 



E 



7r x 7T „ ' — ' nix — in) ' 

0^n=— 00 v ' 



(3.44) 



one finds the poles 



At n = 4 ir n i Rh , 



(3.45) 
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for any integer value of n. For Q > (corresponding to emission) the contour of integration can 
be closed along an arc in the lower half plane and we get 



P(2fi;ti,L) ~ N 2 ^^ [ L dT \ A ^ y e -an*nR Hy /rHtejR 

V ' 2m 2 Jo R?Wy/\-R H /R ^ 

+P L , (3.46) 

where Nl is the maximum value of n for which the pole At n lies within the contour of integration 
and Pl is the contribution from the arc. 

If Ni is sufficiently large, then Pi can be neglected, 

E e-^^V^m „ * , (3.47) 

n=l 

and one recovers the Planck distribution with (instantaneous) temperature 

T = - Th : . (3.48) 
y/1 - Rh/R 

In general, however, one expects that the time L is not long enough and non-thermal contri- 
butions will render the evaluation of the probability of emission very complicated. One should 
then divide the time of collapse into small intervals and compute all the relevant quantities step 
by step. This procedure can be simplified by employing an approximation introduced in order 



to follow a trajectory [35] and which amounts to estimating the probability of emission per unit 

(3.49) 



time at a given value of t (or r) as 

dP(Q;t) Tre 2 h A^ 2 |A ^ 2 | 2 



dL 2m 2 RSWy/l- Rh/R e 2 h n / k B t - 1 ' 

in which all time-dependent quantities must be evaluated at the same time t and the right hand 
side is to lowest order in R. 

We finally recall that in Ref. || we proposed that the finite thickness of the macroshell could 
provide a solution to the problem of ultra-Planckian frequencies. In fact, if 10* is the frequency 
of the emitted quanta as is measured by a distant observer, a fixed observer located near the 
point of emission at r = R will instead measure a blue-shifted frequency 

^=(l-^f) V2 , (3.50) 

and this expression clearly gives u > l/£ p for R sufficiently close to Rh- In order to probe these 
modes, one should use a detector localized in a region smaller than uj~ 1 ~ £ p , or, alternatively, 
one expects that the spread 5 of the wavefunctions $ ra should be less than £ p for our collapsing 
microshells to couple with conformal quanta of ultra-Planckian energies. However, as we have 
shown in Refs. || ; taking 5/£ p — > decouples the emitter from the radiation and ultra- 
Planckian frequencies are not excited (for a similar argument against ultra-Planckian frequencies, 
see Ref. M). 
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3.3 The flux 



It is now straightforward to estimate the total flux of emitted radiation (the luminosity) as 
a function of time. In particular, the rate of proper energy lost by the macroshell per unit 
(Schwarzschild) time of a static observer placed at large r is given by 

— = -A}_^^{uj)T(u)nuj — — — , (3.51) 

where A = Air R 2 is the surface area of the shell, n(u) = (1 - R H /R) 3/2 uJ 2 the phase space 
measure for photons of frequency oj (measured at the shell position r = R), T ~ 1 the grey-body 
factor for zero angular momentum outgoing scalar waves [41| and the last factor is obtained 
from the approximation ( 3.4S| ). The sum in Eq. ( 3.5 1| ) is dominated by the contribution with 
oj = 2 0, and one obtains 



dE 16vr 2 e 2 N 2 £l f S W R H \ R 



in 



(3.52) 



dt 9 R A \Rh J V R ) e 2£lyJl-R H /R/k B T H _ 1 

We can now substitute the expressions ( |3.6| ) for and ( |2.12| ) for 5 and (minus) the flux becomes 
dE 16vr 2 e 2 N 2 ^ R H \ R 2 (3 _ 



dt 9 R 10 ^R A I / 3 \ R) e 2n^i-R H /R/k B T H _ x ' 

We note that, for R close to Rh the Boltzmann exponent is very small (because of the 
Tolman factor). Upon expanding the exponent and keeping the leading order in R — Rh we 
obtain 



**~_f.£(*.Y7l_*?V te. (3.54) 



1/2 

Itt " W \RhJ V ^R 

The flux should therefore vanish for R — > R H and, on comparing with the general relation given 
in Ref. pi, 



dE = Rh ^1-Rh/R + R 2 ~R ^Rh (3 55) 

dt 2G s[l - R H /R + R 2 2G ' 

one obtains, in the limit for R close to Rh, 

Were one to trust the above expression, one should then require that the energy-momentum 
tensor of the radiation remain locally finite when the shell approaches the surface R = R H . 
This amounts to the condition 



Ru - ( 1 - , with 7 > 1 . (3.57) 
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radiation wavelength A 


1(T iU -=- 1(T 7 cm 


S/X 


KT 3 -=- io -1 


R 


10" 3 -f- 10" 2 


total number of emissions 


10 3y 


N s 


10 a& ^_ 1Q w 


total emitted energy 


5% of M s (0) 



Table 1: Typical values of the relevant quantities resulting from the numerical analysis for 
N = 2 ■ 10 40 microshells [R H (0) ~ 4 • 1(T 12 cm]. N s is the average number of emissions while 
the shell radius shrinks a space S. 



Since the quantity B ~ 1 near Rh, Eq. ( |3.55| ) would lead to the conclusion that the velocity of 
the macroshell satisfies 



« 2 ~U-^) > with/3>i (3.58) 



which means that the proper velocity R should vanish at the horizon. Of course, the number of 
approximations we have employed does not allow one to consider the above a rigorous argument, 
but just a suggestive result. In particular, if the emission near the horizon is not thermal [Nl in 
Eq. ( |3.47 ) is small and just a few of the poles are included], the flux satisfies the condition ( |3.57| ) 



for any finite value of R, the backreaction is reduced and our results are significantly modified 
only very close to the horizon where in any case our approximations break down. Indeed, we 
are only able to perform the analysis for R greater than a few times Rh and, in all the cases we 
consider, the microshells are confined with 5 <C Rh- 



3.4 The trajectory 

One can now integrate numerically the equation of motion for the radius of the shell, 

■ o R 2 u Rh G 2 M 2 , . 

together with Eq. ( |3.55| ) for Rh, 

128^ 3 GN 2 fJ? / R H \ R 2 



RH "~^ a R^Rf e 2n^- RHlR/kB T H _ l ' (3 ' 60) 

where a = e 2 /Air. 

Given our construction, the trajectories depend on the parameters N, a and the initial 
condition R(0). In general, one finds that the non-adiabatic excitations occur relatively close to 
the horizon and there is no strong dependence on the initial value of R. In Fig. || we display 
the comparison between the trajectories computed from the above equations, with N = 2 • 10 40 
and various values of a, and that with constant M s = M (see also Table [l]). Owing to our 
approximations, we are not able to display any reliable results as it! gets close to Rh- The 
behaviour of Rh is shown in Fig. [7|. 

Increasing N requires an increase in a in order to preserve the difference between the radi- 
ating and the non-radiating trajectories (see Fig. || for a tenfold increase in N with respect to 
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(a = a ) 



10000 20000 30000 40000 50000 



(a = 4a ) 




(a = 16 ao) 



150000 200000 

r 



200000 400000 600000 800000 



(a = 64 ao) 



Figure 6: Trajectories of the radiating shell R = R(t) in units of Rh{0) with N = 2 • 10 40 and 
four values of the coupling constant a (upper curves) compared to the non-radiating collapse 
(lower curves). All trajectories are evolved from R(0) = 1000i?f/(0). The time r, in all plots, is 
expressed in units of 2 • 1CT 38 iVGeV" 1 with Ti = c = 1 and a = 1.6 ■ 1(T 33 . 
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Rh 



5000 10000 15000 20000 25000 30 



(a = a ) 



Rh 



10000 20000 30000 400UO ^50000 



(a = 4ao) 



Rh 



50000 100000 T3^tr^— 



(a = 16 ao) 



Rh 



200000 400000 600000 S- 



(a = 64 ao) 



Figure 7: Behaviour of the gravitational radius for the radiating shell trajectories in Fig. [||. 
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5000 10000 15000 20000 25000 



200000 300000 



(a = 6.4 • 1(T 31 ) 



(a = 6.4 • 10"^) 



Figure 8: Trajectories of the radiating shell with N = 2 • 10 41 and two values of the coupling 
constant a (upper curves) compared to the non-radiating collapse (lower curves). 
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Figure 9: Velocity of the radiating shell for the fourth case of Fig. |6| and second case of Fig. |S| 
(upper curves) compared to the non-radiating collapse (lower curves). 



Fig. ||). In fact, from Eq. (3.60), on setting R ~ Rh oc N m one finds (apart from dimensional 
constants) 



Rh oc 



a 



AT8/3 m 22/3 ' 



(3.61) 



which gives an estimate of how the effect scales with respect to the number of microshells, their 
mass to and the radiation coupling constant. 

In Fig. H we plot the velocity R of the collapsing shell (together with the velocity of the 
corresponding non-radiating shell) for two relevant cases. It is clear that R remains (negative 
and) small (within a few percent of the speed of light), thus supporting our approximation 
scheme within which we just retain the lowest order in R. In Fig. |l(] the wavelength of emitted 
radiation quanta is plotted as a function of the time in order to show that it remains larger than 
5 (as given in Table |l]). 

Finally, we have also checked that on replacing the Bose-Einstein factor in Eq. ( |3.60| ) with 
the sum of just a few poles [see Eq. (|3.47|) 1 the evolution of the system is not substantially 
changed (only the backreaction is somewhat reduced). 
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Figure 10: Wavelength A of the emitted quanta and ratio A/5 for the fourth case of Fig. ||. 

4 Conclusions 

In this paper we have analyzed the gravitational collapse (in the semiclassical approximation) 
of a macroshell built up of a (large) number of bosonic microshells (s-wave particles). Starting 
from the classical equations of motion, we obtained the potential and an effective Schrodinger 
equation for each microshell. The potential appearing in this equation ensures that the thickness 
of the macroshell does not increase significantly for a long enough piece of the trajectory, thus 
allowing one to take the thin shell limit at some point in the calculations. The time dependence 
of the potential leads to the non-adiabatic excitations (which we computed to leading order in 
the velocity of the macroshell) of the microshell bound states. On coupling the microshells to a 
scalar radiation field, we showed that the increase of kinetic energy can be radiated away and 
the emitted radiation is approximately thermal with a temperature given by the Tolman shifted 
instantaneous Hawking temperature (whether this mechanism can be related to the Hawking 
effect is a point which requires further study). 

The effect is intrinsically quantum mechanical, since it is a consequence of the quantum 
mechanical-bound state nature of the macroshell and the coherence of the emitted radiation. 
This can cause the shell to lose enough energy so that the backreaction on the trajectory of 
the radius is large. For instance, for the case presented, the shell approaches the gravitational 
radius in a time which is more than an order of magnitude longer than the time a non-radiating 
shell would take to cross the horizon, and loses about five percent of its ADM mass as a burst 
of radiation (which is suggestive of the observed gamma-ray bursts |^3j ) . One may also wonder 
at this point whether it is possible for the collapsing matter to produce enough radiation to 
reduce the energy of the shell so as to prevent the formation of an event horizon. Thus, if the 
ADM mass decreased fast enough, the shell itself would never reach the gravitational radius, 
unfortunately our approximations become unreliable as it is approached. 

The model we have proposed clearly hinges on the formation of a condensate. Let us give 
some arguments in favour of this: of relevance for the latter is the fraction of microshells (par- 
ticles) not in the condensate (ground state). This fraction is related to (T/T c ) u , where T c is 
the critical temperature for Bose-Einstein condensation and the (positive) exponent v depends 
on the confining potential and the dimensionality of the system []. For our case of about 10 40 

6 For example, v — 3/2 for a three-dimensional box and v — 3 for a three-dimensional harmonic trap (for a 
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particles and on setting T = Th (Hawking temperature) one has that Th/T c < 10~ 10 which 
renders our description in terms of a condensate plausible. 

Naturally, the relevance of our bosonic model for gravitational collapse which is expected 
to involve ordinary incoherent matter may be questioned. Moreover, it is clear that fermionic 
microshells would lead, because of the Pauli exclusion principle, to a much wider macroshell 
with and a loss of coherence. However, for a mixed fermion-boson system, one may just consider 
the bosonic part and apply our considerations to it since it will condensate in the lowest state. 
Certainly, the presence of fermionic and/or incoherent matter will reduce the effects we have 
illustrated. Nonetheless, one might even argue that an analogous phenomenon can happen for 
all collapsing bosonic matter, including the bosonic fraction of accretion disks around black 
holes. Our result would thus suggest a new mechanism by which the accreting matter can emit 
radiation. 
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A Microshell distribution and effective potential in NR 

The NR for each microshell of rest mass m starting at large radius (fi 3> 2G M s ) with negligible 
initial velocity can be obtained by setting M> — M < = m and keeping terms up to first order in 



G and to leading order in m in the effective Hamiltonian constraint (2.4). The total energy of the 
system is thus given by the sum of the kinetic energies of all microshells and mutual interaction 
potential energies. The latter can be written, on introducing relative (fj) and centre-mass (R) 
coordinates, as || 

M 2 m 2 N 

V(n, ...,r N ) = -G— + G^J2\fi- rj\ ■ (A.1) 

i<j 

Unfortunately, this potential does not allow for an analytic treatment of the complete problem, 
therefore we shall apply an analogous procedure to that followed in Section ^ in order to obtain 
a Schrodinger equation for each microshell in the mean field of the others. 

Let us first try with a configuration {rj,Mj}/j 0m in which N — 1 microshells have equally 
spaced radii, rj = r\ + (i — 1) 8/(N — 1), so that Mj = m [(r* — ri) (N — l)/8 + 1]. Hence, the 
equation of motion ( |2.3| ) in NR for the N th microshell with radius r\ < x < rjv-i becomes 

1 / dx\ 2 „m 2 — r\ A , . „. 

2 m u) - G vr— +i ) • (A - 2) 

where t = t x is the Newtonian time. Changing to the relative coordinate x = x — R, one obtains 

1 / dx\ 2 ^mM x f x 



H ^2 m {n) +G — Ar- 1 ) ' (A ' 3) 

where we have used |x| <C R and omitted all potential terms which do not contain x (they 
would just contribute to the ground state energy in the Schrodinger equation). It is clear that 

review, see Ref. @). 
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the potential above for x = ±6/2 does not match smoothly with the NR limit of the external 
potential given in Eq. ( |2.9| ) because of the term linear in x. This means that the configuration 
{rj, Mi} is not stable. In fact, the minimum of the potential is not at the origin x = 0, and 
one then expects the microshells to move away from the above configuration. 

In order to obtain an internal configuration compatible with the external potential, let us 
assume a potential with a minimum at x = (as is required by the symmetry of the external 
potential) and consider small oscillations about this equilibrium point. Such oscillations will be 
harmonic and the potential to lowest order is given by 



V NR = G^- § (x-R) 2 + C 1 



(A.4) 



where C\ is a constant. If we assume Mj+i = Mj + m, for the most likely configuration {n, Mi}, 
Ti can be determined by comparing the potential Vnr with the one obtained from the NR limit 
of Eq. (P^) , to wit 



Vnr 



c m(m-2Mj) | 
2 x 



n < x < n+i 



(A.5) 



where C2 is another constant. On equating the potentials Vnr and V^ R in the limit x — > 
(from above) and setting C = C% — C\, one obtains 



2 Mi 



m 



2C 



M 



Gm R 2 5 



R 



(A.6) 



Then, on defining r^+i = r j + yi and subtracting the same expression for x 



r i t fl , one has 



2 m 



2M i+1 
2C 



Gm 
2C 
Gm 



y% + 



2 Mi 
Mrf 



JJi 



(r+-R) 2 -(r+ + 



R 



ViM 
R 2 5 



(• 



R 2 5 



+ Vi-R 



R 2 -4r+R + 3 (r+) 2 + Vi (Zrf -2R + yi 



(A.7) 



The first term on the right hand side corresponds to the homogeneous distribution if the con- 
stant C = Gm 2 (N — l)/6. In general, however, the other terms cannot be neglected, so 
that the microshells are not equally spaced. Further, the yi can be computed recursively for 
i = 1, 2, . . . , N — 2 on making use of Eq. ( |A.7|) and the condition t\ = R — 5/2 (which deter- 
mines the constant C). In any case, as we discuss in Appendix]^, the actual quadratic or linear 
structure of the potential inside the macroshell has little influence on the energy levels of the 
lowest states. 



B Shell dynamics in Schwarzschild time 

It is known that, even if two theories are classically equivalent, their quantization may lead 
to different Hilbert spaces and, therefore, different quantum pictures and that this problem is 
potentially present in Einstein's gravity (see, e.g., Ref. |30f| ) . It is thus interesting to compare the 
potential Vm^ obtained in Section [2] with the one corresponding to a foliation of space-time into 
slices parameterized by the Schwarzschild time = tjy measured by a (distant) static observer. 
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Given this potential, one can (in principle) construct the Hilbert space 7i T and compute any 
observable in order to check the equivalence of H T and Hf Unfortunately, the form of the 
potential we obtain in the following makes the explicit construction practically impossible. 

First we observe that the proper time of the test shell can be expressed in terms of the inner 
and outer Schwarzschild times t < and i> = t x (see Fig. ffij) by making use of 



dri 



which yields the ratio 

dt< 

where, from Eq. (EO), 



2 G M. 



x 

2GM~ 



x 



dti 



dti 



2GM. 



x 

2GM 



-i 



> 



x 



dx' 



dx' 



x-2GM > 2x 2 G(M > - M<) f dx 

x-2GM < ~ (x -2GM > ){x -2GM < ) 2 \dt^ 



dx 
dtZ 



1 (x-2GM > ) 2 F(x) 
x x-2GM > + xF(x) 



(B.l) 



(B.2) 



(B.3) 



These expressions can be iterated for each microshell with radius x < < and, together 
with (GO), yield 



dr x 

dtnr 



1 



2GM X 



N 

>< n 

i=x+l 



2GM;, 



2 G M T 



dx\ 2 
~dTx) 



2r 2 i G{M i -M^ 1 ] 



n-2G (r t -2GM i )(r i -2G 



<h_ ^ 

dt, 



1 



2GM X 



+ 



(x-2GM x )F(x) 
x - 2 G M x + x F(x) 



N 

>< n 

i=x+l 



2G Mi 2rfG(M i -M i _ 1 )(r l -2GM l )F(r l 



n — 2G Mi-i (n - 2GM W ) 2 [n-2GMi + nFin)} 



(B.4) 



where r x+ \ is the radius of the first microshell having radius greater than that of the test shell 
and 



-1 + 



{M t -Mi - ^ 



+ G 



Mi + Mi_i 



m 



+ G 



(B.5) 



On then multiplying Eq. (|2.3|) by m/2 and using Eq. ( [B.4D , one obtains the equation of motion 
for the microshell in the form of an effective Hamiltonian constraint in the Schwarzschild time 



dx 
dr x 



dr x 
dtoc 



+ 



V® = . 



(B.6) 



The explicit form of 7" now depends both on the distribution {r^, Mi} of the microshells inside 
the macroshell and the cumbersome time conversion factor in Eq. (B.4p. 
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Hence, we again rely on the thin shell limit (<5 <C ri) and consider the form of the potential 
when x < r\ or x > rjv- For x < X one finds [see case (1) in Fig. ^] 



dx 

dtnc 



dX 



2Gm {x-2Gm)F < (x) 



+ 



x x — 2Gm + xF < (x)_ 
X-2GM S 2XG(M s -m)(X -2GM S )F < (X) 

_ X -2Gm + (X -2Gm) 2 [X - 2GM S + X F K (X)] 



F<(x) 



(B.7) 



2GM S (X - 2GM S )F < {X) 



+ 



X X-2GM S + XF < {X) 



F<(X) . 



Analogously, for x > X [see case (2) in Fig. |2| 



dX 



2G(M s -m) [X -2G(M s -m)]F > (X) 



+ 



X ' X - 2G{M S -m) +XF > {X) 
x-2GM s 2xGm(x - 2GM s )F > (x) 

_x-2G(M s -m) + [x-2G{M s - m)] 2 [x - 2G M s + x F > {x)} 
x F>(X) 



(B.8) 



dx 

dtrsc 



i 2GM S + (x-2GM s )F > (x) 



x -2GM S +xF > (x) 



F>(x) 



where F < and F > have been given in Eqs. (2.5) and (|2.£|). On introducing the relative (f) and 
centre-mass (R) radii, after some lengthy algebra one obtains an effective Hamiltonian for the 
two shells given by 



2 \dtoo. 



1 



dr 



2 r ~ \dt c 

where the potential can be expanded in powers of G and f/R according to 

OO OO > _ v jfe 

V (t) = G n Y, V n,k(M,rn;M s ;R) [-) . 



n=0 



k=0 



(Bi 



(B.10) 



Only the terms linear in r/R are relevant for our calculation and, to leading order in G, they 
are given by 



M s m 
2R 
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, M 2 M? M 4 



5 

f>+- 
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(B.ll) 



where we have again neglected non- leading terms in m. 

On interpolating between the two terms given by V\,\ in Eq. we obtain 



(B.12) 
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with the potential 



2R 2 d 



2-rr + l + 3 — 



M M r 



M 



2ii 2 



/ M 2 M s M M 2 M 3 M 4 \ 

^- 2 M- 18 + 3 MT + 6 M?"Mf"M? ' (R13) 



At) 



which holds for \r\ < 5/2. 

If we now take the thin shell limit of Vm only the linear part V\ \ survives. We can then 
compare the result with the corresponding linear part of Vm \ as computed in Section |], using 
the Schwarzschild time, and then obtaining 



v {r) (t) = y (r 



dr 



where 



and 



dr 

dtnc. 



( 2GM S + (R-2GM S )F{R) 



R 



R-2GM S + RF(R) 



F(R)- 

This yields, to lowest order in G, 



M 2 M s G 2 M 2 

1 + W + G li + ^W 



Vi T) (t)^G 



M s m \r\ 
2R 2 



M 2 
M 2 



(B.14) 



(B.15) 



(B.16) 



(B.17) 



which coincides with the thin shell limit (i.e., linear part) of Vm \ Eq. ( |B,13; ), for M = M s . This 
shows that the effective potential governing the motion of the microshells is the same in NR for 
the two foliations {S^} and {S T }. 



C Microshell confinement 

We shall now show that the expression for the macroshell thickness 5 in Eq. ( |2.12| ) can be 
obtained within the adiabatic approximation (R constant) in NR by either considering just the 
linear part of the potential Vm ^ in Hm (thus neglecting the details of the distribution of the 
microshells) or just the quadratic part (which is related to a given distribution of microshells, 
see Appendix [A]) . From the Hamiltonian H$ one can obtain a (time- independent) Schrodinger 
equation 



7T" 

2 m 



+ V ( 



(C.l) 



where tt t = —ifidf and, on retaining the linear part, one has (N m = M = M s ) 



Vi 
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GNm 2 
2R 
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Figure 11: The potential Vjh V at order G 2 for M = 10 40 m h and 5 = i h = 1CT 16 m 1 < M/M s < 2. 
Vertical units are arbitrary. 



or, on considering the quadratic term, 



(2) 



GNm 2 f 2 
2R 2 J 



(C.3) 



Solutions to Eq. ( C.l ) with Va-\ = Vm are given (apart from a normalization factor) by the Airy 
function Ai [45|, 



Ai(£-£ 2fe ) 

sgn(f) Ai(£-£ 2 fc+i) 



n = 2k 



n = 2k + l 



(C.4) 



where k E IN" and £ = (2G N m 3 /fi 2 i? 2 ) 1 / 3 f. The £2*1+1 are the zeros of the Airy function, 
Ai(— ^2fc+i) = 0, and ^2k the zeros of its first derivative, Ai'(— £ 2 fc) = 0. The corresponding 
eigenvalues are given by 



£(1) 



N 



1)2 \ 2 / 3 



m 



R 2 



2V3 ' 



(C.5) 



and the spread of such states (for n small) is well approximated by Eq. ( 2.12j ). 

Solutions for VtA = Vn) are given instead by the usual harmonic oscillator wavefunctions 
(again we omit normalization factors), 

^=H n (C)e~^ , (C.6) 

where £ = (G N m 3 /2 fi 2 R 2 <5) 1//4 f and H n are Hermite polynomials [45|. The energy eigenvalues 
are 



! hNm l v 
25 ~R 



n + 



1 
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Figure 12: Plot of (continuous lines) and ^ffl (dotted lines) for n = 0, 1. Vertical units are 
arbitrary and d should be 5. 

If we now substitute the expression for 5 previously obtained, we get the energy spectrum 

which is of the same order of magnitude as Eq. ( |C.5| ). Also, the spread of such for n small 
is again given by Eq. ( |2l^ ) (see also Fig. 0) and we can therefore represent the state of each 



microshell either by using 

It is also important to observe that the energy gap between two adjacent states, 

AE ~[-W) m ~l-R^J m ' (C - 9) 

is much smaller than m for R > Rh, which suggests that outside the gravitational radius the 
effects considered in the present paper dominate over the microshell creation studied in Ref. §. 

We may also see the small difference between the energy levels in the two cases by employing 
a W.K.B. quantization K|. Indeed if we require 

7r r df = m <j> f df = + — ^ h , (C.10) 

for both potentials and consider cycles within the potential well, that is the amplitude of the 
oscillations will correspond to the width of the shell, one has 



hS/2 

rn<t>fdf = V2ml df 

8/2 



V {l) (5/2)-V {t) (r)] 1/2 = {n+ l ^h, (C.ll) 
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thus determining 5^ and <W f° r the linear and quadratic potentials respectively. One obtains 



5 (D 
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2 ( n+ 2) k 16GMra 



1/2' 
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and correspondingly a binding energy 
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(C.13) 



which are of the same order of magnitude and we have added to E n the additional 0(G) term 



in Eq. ( 2.1 1| ) so as to have continuity of the potentials at |r| = 8/2. It is clear that in both cases 



we have allowed for a maximal oscillation of the test microshell and thus obtained <5(j) . 

Finally, let us now check how close to the gravitational radius we can use the bound states 
obtained above. From Eq. (|2.9D one finds that is always positive for f < 0, but becomes 
negative for f > if 



R < 



G 2 M 2 
2Rh 



Rc 



(C.14) 



(i) 

As a sufficient condition for applicability of the modes § n , we shall require that the energy be 
less than the value of the potential for f ~ 45 (see Fig. 12), i.e., En < V(i){f = 45). The above 
condition implies n < 4 for R > Rh (> Rc) an d is therefore satisfied for the second excited 
state (n = 2) we use in Section |3|, otherwise the shell breaks up. 
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